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Abstract 


We consider the first pion excitation as a sub-continuum resonance in the 
pseudoscalar channel, and we obtain parameters characterizing this resonance 
through a global fit of the Borel-parameter dependence of the field-theoretical 
pseudoscalar Laplace sum rule to its hadronic (pion -|- pion-excitation -|- QCD- 
continuum) content. Our analysis incorporates finite-width deviations from 
the narrow resonance approximation, instanton effects, and higher-loop per¬ 
turbative contributions to the pseudoscalar correlator. We obtain the follow¬ 
ing values (uncertainties reflect 90% confidence levels): mass Mn = 1-15 ± 
0.28 GeV, width Pn = < 0.48 GeV, decay constant r = [FnM^ / firm^]^ = 
4.7 ± 2.8. 
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The mass, width, and decay constant of the first pion excitation are of interest to effective 
theories of low energy QCD. Theoretical argnments have been advanced to motivate masses 
below 1 GeV for this state,despite the Mn = 1300 ± 100 MeV entry in the present 
edition of the Particle Data Gnide.^ The possibility that this resonance may be very broad 
is also not excluded by present data, which gives a width Tn between 200 and 600 MeV.^ A 
broad width is suggestive of a small value for this resonance’s decay constant Fn, which is 
experimentally undetermined. Recent theoretical work, however, suggests that Fn may be 
very small, even less than 1 MeV, if Mn is as large as 1300 MeV.^ If the first pion excitation is 
as light as 1 GeV, one cannot automatically assume that this resonance’s contribution to the 
correlation function of pseudoscalar currents is absorbed in the continuum of hadronic states 
approximated by purely-perturbative QGD,^ a continuum whose onset is generally assumed 
to be substantially above s = 1 GeV^. ^ Gonsequently, a QGD sum-rule approach® to the 
1=1 pseudoscalar resonances cannot necessarily be limited to the lowest-lying resonance 
(the pion). 

The idea of applying QGD sum-rule methods to the properties of the first pion excitation 
is an old one. Hubschmid and Mallik^ included the hrst pion-excitation state within a hnite- 
energy sum-rule analysis of the pseudoscalar channel, assuming a mass Mn = 1.5 GeV 
based upon a linear Regge trajectory with a universal slope. Their results are suggestive 
(up to large uncertainties) of a substantially larger value for Fn than anticipated from 
phenomenological hadronic physics. Their results, however, are obtained within the context 
of a narrow resonance approximation in which Fn is assumed to be zero, as their analysis was 
directed more towards a determination of quark masses than of parameters characterizing 
the hrst pion excitation. 

The work presented here extracts information about the hrst pion excitation’s mass, 
decay constant, and width from a global ht of the Borel transform of the pseudoscalar- 
current correlator to its hadronic content, assuming that both the pion and the first pion 
excitation reside below the continuum threshold. Subsequent pion excitations in our anal¬ 
ysis are assumed either to be above the continuum threshold, or to have decay constants 
sufficiently small to decouple them from sub- continuum physics. We depart from prior anal¬ 
yses 6F,8,9, 10,11 pseudoscalar channel sum rules by incorporating 1) direct instanton 

contributions, 2) hnite-width ehects, and 3) higher-order perturbative contributions to the 
pseudoscalar current correlation function: 

1. Direct Instanton Contributions: In the QGD sum-rule approach, long distance ehects 
are characterized by local matrix elements of quark and gluon operators averaged over 
the physical vacuum - the QGD-vacuum condensates. However, such local condensates, 
corresponding to vacuum huctuations of inhnite correlation length, are insufficient to account 
for the full nonperturbative content of the QGD vacuum, as they do not take into account 
vacuum huctuations of non-local origin arising from instantons.^^ Such hnite-correlation- 
length huctuations may provide the underlying mechanism for chiral symmetry breaking, 
and need to be separately included in an analysis of either scalar or pseudoscalar current 
correlation functions.^® The direct single-instanton contribution to the Borel transform of 
the pseudoscalar correlation function has been calculated in the instanton liquid model: 
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with instanton size denoted by pc, and the light quark mass in the SU(2)- symmetry limit 
denoted by m^. Note that instanton effects can be “turned off” by letting oo. In eq. 
(1), also note that n(r) is the Borel transform to the 1 = 1 component of the pseudoscalar- 
current correlation function: 

n(r) = - r Im[UP{s)]e-^^ds (2) 
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2. Finite- Width Effects: Conventional methodology for determining the hadronic content 
of QCD sum-rules rests upon the narrow resonance approximation, in which hadronic Breit- 
Wigner contributions to the imaginary parts of appropriate correlators are proportional to 
(5-functions:® 
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One can undo the narrow resonance approximation through explicit utilization of the Breit 
Wigner shaped^ Analysis of the longitudinal component of the axial-vector current correla¬ 
tor [corresponding to dm^/s times the pseudoscalar current correlator] yields the following 
hadronic contribution: 
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The parameter sq represents the continuum threshold above which hadronic physics coincides 
with purely-perturbative QCD. In the present work, we hrst express the imaginary part of 
the Breit-Wigner shape occurring in (6) as a Riemann sum of unit-area rectangular pulses 
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A semi-analytic approximation for the integral over the Breit- Wigner shape can be realized 
by noting (for sq > -|- MnT') that 
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In the present analysis, we have utilized a 4-pulse approximation to the right-hand side of 
(8), choosing f = 0.7 in order that the area of the four pulses be equivalent to the total 
area (= tt) under the Breit-Wigner curve. [In the n ^ oo limit, such equivalence is true for 
any choice of f between 0 and 1]. We then obtain a simple expression for the hnite-width 
correction W to the narrow-resonance-approximation contribution: 






ds = 


"n-W[Mn,rn,r] 


W[M, r, r] = 0.5589A(M, 3.5119 B, r) + 0.2294A(M, 1.4412 B, r) 

+ 0.1368A(M, 0.8597 B, r) + 0.0733A(M, 0.4606 B, r). (11) 

Note that hB —1 in the narrow-resonance limit (Bn —^ 0). The four- pulse approximation 
we utilize serves to mitigate numerical difficulties that arise from the inhnite Breit-Wigner 
tail, which extends unphysically into both Euclidean (s < 0) and post-continuum (s > Sq) 
values of s. 

3. Higher-Order Perturbative Contributions: Virtually all previous sum- rule treatments 
of the pseudoscalar channel have made use only of the one-loop expression [or incomplete 
2-loop expressions] for the purely perturbative QCD-contribution to the pseudoscalar cur¬ 
rent correlator. However, the perturbative content of the pseudoscalar current correlation 
function has been known for some time to three-loop order, with higher-loop contributions 
accompanied by numerically large coefficients:^^ 

Im[UP{s)r^^ = ^ - 21n(sr) 

87r I TT L 3 

+ ^ [40.684 - 31.667In(sr) - lA177r‘^ + 4.251 ln2(sr)] + 0{al) | . (12) 

The renormalization point /i = is chosen to identify the Borel parameter as the 

renormalization-group (RG) scale parameter. In the work presented here, we have con¬ 
sidered both 2- and 3-loop expressions for n^(s) in determining parameters characterizing 
the hrst pion excitation. Computational limitations prevented us from determining 90% 
conhdence-level bounds on our htted parameters when using the 3- loop expression-we 
could obtain numerically efhcient Monte Carlo estimates of these conhdence levels only for 
the 2 loop case. We were, however, able to compare explicitly the y^-minimizing hts ob¬ 
tained when either 2- or 3- loop expressions for n^(s) are utilized. Rows 5 and 6 of Table I 
demonstrate that the pion-resonance mass, decay constant, and width are shifted very little 
in going from two to three loops, despite potentially large 3-loop perturbative corrections. 
It should be noted that the minimum increases by an order of magnitude when the pion 
resonance is omitted (leaving only the ground state pion). This indicates that, unlike the 
case for other channels, the pion resonance is too strong to be absorbed into the continuum. 
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Our analysis of pion-resonance properties follows from a weighted least squares fit of 
the r-dependence [0.4 GeV“^ < t < 2.5 GeV“^] of the hadronic contributions to Ri(t) 
to the corresponding field theoretical expression, as derived from the pseudoscalar current 
correlation function. 

The region of r is chosen by placing a 20% npper bonnd on the theoretical uncertainty 
in based npon a 30% continunm and a 50% power-law-nncertainty. The hadronic 

contribntion, which explicitly contains the htted pion-resonance parameters Mn, Fn, Tn and 
the continnnm threshold Sq, is obtained via snbstitntion of eqs (10-12) into (6). The held- 
theoretical contribntion i?(*(r)[= as dehned by eq. (2)] contains direct instanton 

contribntions [eq. (1)], QGD-vacnnm condensate contribntions,®’®’^® and pnrely-perturbative 
contribntions that can be extracted by snbstitntion of (12) into (2). Since the qnantity i?i*(r) 
satisfies a renormalization gronp eqnation,^’^ conpling constants as and qnark masses niq 
appearing ontside RG-invariant condensates can be replaced with scale-parameter dependent 
RG-improved expressions referenced to the /i = mass scale to the appropriate loop 

level. The two loop expressions are: 
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The parameter m is the RG-invariant qnark mass, which we inclnde as one of our htted 
parameters by comparison of the explicit r- dependence of 
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SU(2) breaking effects are higher order in m^, since their contribntions to Ri are in terms 
proportional to {mu—mdYirriu+mdY-^ The contribntions to Ri from the dimension-5 qnark- 
antiqnark-glnon condensate and the dimension-6 triple-glnon condensate are proportional 
to and are therefore not inclnded in (17). Onr £t is obtained using a standard set 

of pertnrbative and nonpertnrbative RG-(quasi-)invariant QGD parameters: < mqq >= 
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-f^ml/A [U = 131 MeV], < > = 0.045 GeV^ A = 0.15 GeV, pc = (600 MeV)-^ Input 

assumptions concerning the deviation of the dimension-6 fermion condensate < Ois{Q<iY > 
from its standard® value 1.8 ■ 10“^ GeV® = Sq are specihed in Table I, and are discussed 
below. Optimal values for the parameters a, r, Mn, Tn, and Sq correspond to the £t of the 
r dependence of (15) to (17) that minimizes a weighted for the previously described 50% 
uncertainty for power-law corrections and a 30% uncertainty for continuum corrections. As 
mentioned earlier, uncertainties in the fitted parameters are obtained only for the two-loop 
case [for which the integral in (15) can be easily determined analytically] from a Monte- 
Garlo simulation incorporating the power-law and continuum uncertainties described above, 
as well as a 15% variation in the value of pc, and (in rows 1-4 of Table I), a factor of two 
vacuum-saturation uncertainty in the value Sq for < as^qq)"^ >. 

The value Sq was originally obtained ® through use of vacuum saturation [< (gg)^ >=< 
gg >^] in conjunction with the value for < mqq > appropriate for a 5 MeV quark mass. 
Either a heavier quark mass and/or a direct violation of the vacuum-saturation assumption 
contribute substantial uncertainty to this quantity, which our Monte-Garlo simulation allows 
to vary over a factor of 2. To demonstrate the overall insensitivity of our fitting procedure to 
the value of < Q;s(gg)^ >, we have constrained this condensate to equal exactly the standard 
value Sq in rows 5 and 6 of Table I. Gomparison of rows 4 and 5 demonstrates that removal 
of the Monte-Garlo variation in the magnitude of the dimension-6 condensate has virtually 
no effect on our results. The results in rows 5 and 6, however, are indicative of a substantial 
violation of vacuum-saturation: the standard parameter value is much larger than that 
anticipated for < as{qqY > via vacuum saturation if the quark mass at /i = 1 GeV is of 
order 8-9 MeV. Row 7 of the Table corresponds to a fit in which exact vacuum-saturation is 
imposed on the value of < ctsiqqY >, by including the implicit dependence of this quantity 
on the fitted RG-invariant quark mass rh. Thus, the requirements < (gg)^ >=< gg and 
< gg >= —f^ml/img are simultaneously upheld. When included in our fitting procedure, 
these input assumptions lead to a p = 1 GeV quark mass much larger than 5 MeV, although 
the upper end of the range obtained for Mn remains consistent with the 1200 MeV empirical 
lower bound.^ 

To conclude, 90% confidence-level results for the first pion-excitation are summarized in 
rows 4-7 of Table I. These results incorporate 1) the explicit direct-instanton contribution, 
2) explicit hnite-width effects, 3) the full two- and three-loop perturbative contribution to 
the sum rule, and 4) the violation [or explicit preservation (row 7)] of the vacuum saturation 
assumption for the dimension-6 fermionic condensate. We believe the results of row 4 to 
be phenomenologically salient, as the fitted range for Mn is found to be fully inclusive of 
the present 1300 ± 100 MeV empirical range.^ The final three rows suggest a somewhat 
lighter pion-excitation than that indicated in the Particle Data Guide, though a mutually 
consistent value of 1200 MeV is not ruled out. 

Although the results of rows 4-7 vary according to the input assumptions employed, a 
number of common features are evident. The decay constant Fn [using PDG conventions 
consistent with G= 131 MeV] is seen to be at least 2.8 MeV, consistent with the larger 
values of Tn already anticipated by Hubschmid and Mallik.^ Although little sensitivity is 
evident to the fitted value of the width (the upper bounds quoted are 90% confidence levels), 
values larger than 480 MeV for Pn appear to be ruled out. The fits we obtain are found 
(surprisingly) to uphold the narrow resonance aproximation (P = 0). The fitted parameter a 
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allows a determination of the RG-invariant quark mass m, and a subsequent determination 
via (14) of the SU(2)- symmetric quark mass at a /i = 1 GeV momentum scale. This value in 
row 4 is less than values anticipated in the absence of instanton corrections (rows 1 and 3). 
A comparison of rows 3 and 4 also demonstrates that direct instanton contributions reduce 
somewhat the value of Mn, but increase Fn- Finally, we note the relative insensitivity in 
Table I to the value of the continuum threshold sq. The analysis clearly favours a threshold 
above for the onset of equivalence between perturbative and hadronic QGD. The viability 
of larger values for sq suggests that smaller values for the pion-excitation decay constant 
Fn could be most easily accommodated by incorporating additional pion resonances into 
sub-continuum physics. 
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Table I 


INPUTS 

OUTPUTS 

Pe(MeU-i) 

D 

H 

aaifcAiMi 

Mn(GeV) 

r 

Fu (MeV) 

a (GeV^) 

mq[fr=lGeV] 

P (GeV) 


oo 

B 


0-2^6 

1.34 ± 0.16 

8.9 ± 3.2 

4.3 ± 1.3 

0.049 ± 0.016 

9.1 ± 1.5 MeV 

- 

a 

1/600 

B 

B 

0-2^6 

1.0 ± 0.25 

5.4 ± 3.8 

6.0 ± 3.7 

0.063 ± 0.035 

8.0 ± 2.2 MeV 


a 

oo 


B 

0-2^6 

1.31 ± 0.25 

5.7 ± 2.8 

3.6 ± 1.6 

0.048 ± 0.014 

9.2 ± 1.3 MeV 

<0.50 

a 

1/600 


B 

0-2^6 

1.15 ± 0.28 

4.7 ± 2.8 

4.2 ± 2.4 

0.069 ± 0.037 

7.7 ± 2.2 MeV 

<0.48 

a 

1/600 


B 

^6 

1.07 ± 0.17 

4.6 ± 1.8 

4.8 ± 1.8 

0.052 ± 0.014 

8.9 ± 1.2 MeV 

<0.48 


1/600 


B 


1.05 

4.36 

4.86 

0.0622 

8.10 MeV 

0 

a 

1/600 


B 


0.95 ± 0.23 

6.3 ± 3.8 

7.2 ±4.1 

0.036 ± 0.019 

10.6 ± 2.8 MeV 

< 0.30 

a 


Fit of n(1300) pion-resonance parameters to 90% confidence levels under various input 
assumptions: 

Pc = oo\ No direct-instanton contribution, 
r = 0: Narrow resonance approximation. 

I = 2: Full two-loop perturbative contribution. 

I = 3: Full three-loop perturbative contribution to the sum rule. For this case, only the £t 
which minimizes is quoted. Computational limitations precluded the determination 
of 90% confidence levels. 

Sq: Sq = 1.8 ■ 10“^ GeV® (see text), corresponding to exact vacuum saturation when rh = 5 
MeV. Row 7 corresponds to exact vacuum-saturation of < as{qqY > for the fitted values 

of rh, TTlq. 


9 





































































